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Abstract 

We present the full NLO QCD corrections to the production of a W^W~ pair in association with two jets 
in hadronic collisions, which is an important background to New Physics and Higgs boson searches. We 
include leptonic decays of the W-bosons with full spin correlations. We find NLO corrections of the order of 
10% for standard cuts at the LHC. The scale dependence is considerably reduced by the NLO corrections. 
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1. Introduction 

The pair production of electroweak VF-bosons 
in association with jets is an important source of 
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backgrounds for new physics searches at the LHC. 
Events stemming from new physics interactions, 
like cascade decays in supersymmetric extensions 
of the Standard Model, will produce energetic lep- 
tons and multiple jets in combination with missing 
energy. The same signature is produced by Stan- 
dard Model (SM) processes involving the produc- 
tion of leptonically decaying M^-bosons and jets, 
where neutrinos escape undetected. A similar sig- 
nature also emerges in top-quark pair production 



within the semileptonic decay mode, and in the vec- 
tor boson fusion (VBF) Higgs-search channels. In 
the latter case, in addition to W-pairs originating 
from Higgs decay, the identification of the process 
relies on the tagging of two forward jets. Therefore, 
the direct production of + 2 jet final states 

constitutes an irreducible background to prominent 
signals. Further, they can serve as a test of the 
electroweak (EW) interactions, allowing to measure 
Higgs boson couplings as well as to probe the pres- 
ence of anomalous vector boson couplings [J, Q . 

Leading order estimates of the background leave 
us with a large uncertainty about the overall nor- 
malisation of the event rates, hence calculations at 
next-to-leading order (NLO) in the strong coupling 
constant are necessary to obtain reliable predic- 
tions of the background from the SM production 
of plus jets. 

NLO calculations of M^-boson pairs in association 
with zero or one jet showed significant higher order 
corrections. The QCD corrections to pp — )■ W^W~ 
have been calculated in Refs. [1, |3l while the EW 
corrections have been computed in [5| . Phenomeno- 
logical studies can be found e.g. in 0, 0, S 0, 
NLO QCD corrections to W+W~ production 
within the Randall-Sundrum model have been cal- 
culated in pp — ^ W^W' plus one jet at NLO, 
including decays to leptons, has been studied in 
Refs. 12, [13, [3li where a K-factor of about 1.3 
has been found. Further, the loop induced process 
gg — ?> W^W^ including leptonic decays has been 
calculated in Refs. 1^ 3 121: for phenomenologi- 
cal studies with emphasis on Higgs search see also 



lOj . It was found that the application of real- 
istic Higgs search selection cuts increases the cor- 
rection from this channel up to 30% of the total 
pp — ?> W^W~ — ?> leptons cross section. 

A VF-boson pair in association with two jets 
can be produced through both vector boson fusion 
(VBF) or production mechanisms without weak 
bosons in the t-channel. NLO QCD corrections 
to the VBF mechanism have been calculated in 
[H, [20|. More recently, NLO QCD calculations 
have been presented for the non-VBF case: The 
production of equal-charge W^-bosons, W^W^jj, 
has been computed at NLO accuracy in 12111. It also 
has been combined with a parton shower [22] within 
the POWHEG framework [23, 24], and later has been 
complemented by EW corrections :25| . First results 
for NLO QCD corrections to the process W^W~jj 
at hadron colliders have been presented in see 
also f27l. 



In this letter we present an independent calcula- 
tion of pp — W~^W~jj, which takes into account 
contributions previously omitted. In particular, we 
also consider diagrams where the W~^W~ pair or 
the Z/j bosons couple directly to a virtual quark 
loop and the contributions from the third quark 
generation in closed fermion loops, and quantify 
their effects on the total cross section. We also 
include the decay of the W bosons into leptons re- 
taining full spin correlation exploiting the narrow 
width approximation. 

In the past years, very rapid progress in the calcu- 
lation of higher order corrections to multi-particle 
processes has been achieved, both due to the emer- 
gence of unitarity based methods and the further 
development of techniques using a tensor reduction 
of Feynman diagrams. As a consequence, NLO re- 
sults for hadron collider processes with up to five 
particles in the final state have been produced in 
the last few years, for processes involving one-loop 
functions with six or more external legs see e.g. 
pp W/Z + 4jets_0, (i?!, pp -)■ 4jcts [30^, pp 
mz^l + ajets (siTia, M, JS, 36],_pp ttbb 
Mm,^,PP^tt + 2 jets Jii , pp bbbb [4l|, \^ , 
pp ti W+W-bb "44 pp ^ W+W+ + 2 
jets [2l|, pp W+W- + 2 ids 26], VK77-l-jet [H, 
e+e" ^> 5 jets [Uli^]. 



The evaluation of the NLO QCD corrections 
for W^-pair production in association with two jets 
can also be considered an important benchmark in 
view of the advances in the development of pack- 
ages for fully automated NLO calculations 48, 5^ 
[Hol , [51I , [52I . As described in the next section, for 
our calculation we employ the combined use of 
MadGraph 53L MadDipole [H, [Hi], MadEvent 
and GoSam [52j. Our results demonstrate the abil- 
ity of the GoSam package to handle the evaluation 
of highly non-trivial one-loop amplitudes in an effi- 
cient and automated way. In particular, the possi- 
bility to filter certain classes of diagrams and pro- 
cess them separately allows us to quantify the rela- 
tive impact of different contributions. For example, 
we find that the contribution from quarks of the 3rd 
generation in the loops amounts to about 3% of the 
virtual contributions restricted to two generations. 
Overall, the NLO corrections at = 2 Mw increase 
the cross section by about 10%, and the scale de- 
pendence is considerably reduced. 

This paper is organized as follows. In Section 
[5] we describe technical details of the calculation, 
while Section [3] contains some phenomenological 
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studies. Finally, in Appendix A we provide numeri- 
cal results for the virtual amplitudes of the different 
partonic sub-processes for comparison. 

2. Technical Details 

2.1. Description of the Calculation 

The leading order (LO) and the real radiation 
matrix elements are generated using MadGraph [ssl 
SgI . 57 1 . For the subtraction of the infrared singu- 
larities we use Catani-Seymour dipoles [H^, sup- 
plemented with a slicing parameter a as proposed 
in [H^ implemented in the package MadDipole [s^, 

13. 

The code for the evaluation of the virtual cor- 
rections is generated by the program package 



GoSam [52| which combines the automated ap- 



proach of diagram generation and processing 
with integrand- reduction techniques 62 , 63l 

'MM 



6C, 
64, 



61 

65 . 66l 51, 63] supplemented by the automatic al- 
gebraic computation of the full rational part. 

For the tree-level, the real emission part and the 
integrated subtraction terms, the integration over 
the phase space is carried out using an in house ver- 
sion of MadEvent ^56|]. We have used 2 • 10^ phase 
space points for the tree-level integration and the 
integrated subtraction terms, and 3-10^ for the real 
emission part. For the the virtual contributions we 
have produced several samples of unweighted events 
according to the tree-level processes and performed 
a reweighting of the events to obtain the virtual cor- 
rections. This has been done independently for each 
subprocess. In total, which means summed over 
all subprocesses, we have used 4 • 10^ unweighted 
events. In order to parallelize as much as possible 
the reweighting has been done separately for each 
non- vanishing helicity. The unweighted events also 
have been produced with MadEvent. 

2.2. Partonic Subprocesses 

The description of LO and NLO corrections 
to pp — )■ W^W~ jj requires two different par- 
tonic structures, namely processes with either four 
quarks or two quarks and two gluons, in addition 
to the ly-pair. 

The cross section for the associated production 
of a ly-boson pair accompanied by 6-jets is domi- 
nated by ti on-shell production and decays. These 
processes have been studied in 0, [3, [il] . Here 
we assume the use of an efficient anti-6-tagging 
to exclude those processes and concentrate on the 



"genuine" , i.e. not top-mediated production of the 
two W-bosons and two light jets. In this sense 
we work in the four flavour scheme and treat the 
bottom quarks as massive particles. We use four 
flavour pdfs and no mixing among the generations. 
The bottom and top quark contributions to the 
gluon self-energy are subtracted at zero momen- 
tum transfer and the running of as is given by 
the four light flavours only. Furthermore, we only 
take into account doubly-resonant diagrams, i.e. 
diagrams with two M^-propagators. The effect of 
singly-resonant diagrams, with only one M^-boson 
radiated off an intermediate lepton stemming from 
a 7/Z-propagator, is strongly suppressed by phase 
space, as has been already noted in ^3i]. We do not 
include any Higgs exchange diagrams nor interfer- 
ences with electroweak production mechanisms of 
WW jj final states. 

We used the GoSam package to automatically 
produce the amplitudes, and chose dimensional re- 
duction as regularisation scheme. MS renormal- 
isation supplemented by on-shell subtraction for 
the massive quarks is also done automatically by 
GoSam. 

The eight basic partonic processes listed below 
are the building blocks from which all the other 
subprocesses can be obtained via crossings and/or 
change of quark generation. We assume the W- 
bosons decay as W^ — 7> I'e e"*" and W^ -> t'^iM i 
but as we use massless leptons, the final results for 
other lepton flavour/charge assignments are identi- 
cal. Six of these subprocesses do not involve gluons 
in initial or final states, namely 



du 


-> 


csW+W- 


u u 


— >• 


ccW+W- 


dd 




ssW+W~ 


u u 




ddW+ W- 


u u 




uuW+W- 


dd 




ddW+ W- 



while other two subprocesses also involve two glu- 
ons: 



u u 
dd 



ggW+W- 
qqW+W- 



In Appendix A we give numerical results for each 
subprocess, labeling each according to the involved 
partons, i.e. with (d, d, g, g) we mean the subpro- 
cess dd^ g g W^{ve e^) W~{D^ fi"). 
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In the calculation of the virtual corrections for 
each of the subprocesses, the Feynman diagrams 
have been divided into three groups (called A, B, 
and C) which we process separately. Group A 
is defined by neglecting all contributions from the 
quarks of the third generation in the fermion loops 
and neglecting contributions in which vector bosons 
are attached to closed fermion loops, which is the 
approximation used in Ref . [26i] . Group B contains 
only diagrams with vector bosons attached to closed 
fermion loops, and group C contains all contribu- 
tions from the third generation quarks in the loops. 
The last two sets of diagrams, which are found to 
be numerically less significant, have not been con- 
sidered previously in the literature. 

2.3. Numerical Checks 

The leading order and the real emission matrix 
elements, which are generated by Madgraph, have 
been compared with the corresponding matrix ele- 
ments generated by GoSam. The infrared singular- 
ities computed with MadDipole have been checked 
with dedicated routines implemented in GoSam. 

Concerning the virtual amplitudes, we have re- 
produced the results of known processes with a sim- 
ilar structure as the one presented here. We re- 
computed uu^ bb and g g ^ bb W~ 
with on-shell W^-bosons and massless 6-quarks [i^ 
starting with off-shell VF's to recover the on-shell 
results from (multiple-evaluation of) the leptonic 
VF-decay modes. The process uu bbW^ W~ 
was also checked considering the VF's directly as 
on-shell particles. 

In the case of decaying M^-bosons, we also 
checked all results presented in [2l| for pp — >■ 
W+W+ + 2 jets and in [26] for pp ^ W+W^ +2 
jets. Further details about this last comparison are 
provided in the Appendix. 

In addition, we have performed several checks on 
the integrated results. For the tree-level and the 
real emission matrix element, we have compared 
two different ways of obtaining the same cross sec- 
tion. We have treated the M^-bosons as stable par- 
ticles, which leads to a reduction of the phase space 
and hence to an easier integration. We have then 
multiplied this result by the branching fractions to 
leptons and compared it with the full result includ- 
ing the W decays where no cuts on the leptons have 
been applied. This check has been performed for 
the most complicated class of subprocesses involv- 
ing two quarks and two gluons and we have found 



agreement of the two results within the statistical 
uncertainty. 

As a test on the subtraction terms, we have 
checked the independence of the "slicing" param- 
eter a for several classes of subprocesses. 

The procedure of reweighting for the virtual con- 
tributions has been checked by a comparison of the 
total cross section with a result obtained using a 
full phase space integration. As these checks are 
very CPU-intensive, we restricted ourselves to two 
classes of subprocesses, namely the case with four 
different quark flavours and the subprocess with 
two u-type quarks and two gluons, the latter pro- 
viding the largest contribution to the total cross 
section. Both results agree within the statistical 
error, which indicates that the phase space is fully 
resolved by the integration routine. For our final 
results, about one permil of the phase space points 
have been classified as numerically unstable and dis- 
carded. 

3. Phenomenological Results 

In this section we present a selection of phe- 
nomenological results for proton proton collisions 
at the LHC at 7TeV. A complete phenomenologi- 
cal analysis and a more detailed study of all the 
sources of uncertainties is beyond the scope of the 
present paper and therefore left for a future publi- 
cation. 

3.1. Setup 

For all the results and distributions shown in this 
section we have used the parameters listed below. 



Parameters 


Mw = 


-- 80.399 GeV 


Vw 


= 2.085 GeV 


Mz = 


91.188 GeV 


Vz 


= 2.4952 GeV 


Mt = 


171.2 GeV 


Vt -- 


= 0. GeV 


Mb = 


4.7 GeV 


Tb -- 


= 0. GeV 


a{Mz 


) = 1/128.802 


Cw 





The weak mixing angle is calculated from the W 
and Z mass. The strong coupling constant and its 
running are determined by the set of parton distri- 
bution functions. We used the MSTW2008 pdf set 
[69| . where the values for as at leading order and 
next-to-leading order are given by 

ois,Lo{Mz) = 0.13355 , as,NLo{Mz) - 0.1149 , 

and the running is calculated at one loop for the 
tree-level result and at two loops for the next-to- 
leading order parts. As we neglect initial state b- 
quarks, we use the N f — 4 version of the pdf set. 
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For the jet clustering we used an anti-/cT algorithm 
with a cone size of i? = 0.4 provided by the FastJet 
package 



3.2. Cuts 

We show results for the LHC at ^ = 7TeV, 
using a set of standard LHC cuts given by 

PTj >20GeV, |?7j|<3.2, Ai?jj > 0.4 (1) 

for the jets and 

PTj>20GeV, |?7j|<2.4 (2) 

for the charged leptons. In addition we impose a 
cut on the missing transverse energy of 



'T,miss 



> 30 GeV. 



(3) 



Neglecting final state &-quarks is only justified if 
the &-jets can be (anti-)tagged. As ^-tagging only 
works efficiently if the jets are produced centrally, 
we do not allow for jets with a very large rapidity. 

3.3. Results 

In Fig. [T] we show the dependence of the total 
cross section on the variation of the renormalisation 
and factorisation scales. The cross section, within 
the cuts given above, has been evaluated for the 
three values oi fiji = fj,p = Mw,2 Mw, iMw. As 
expected, adding the next-to-leading order correc- 
tions strongly reduces the theoretical uncertainties. 
Choosing the scale /i = 2 M\y as the central scale, 
we find 

'll3%{sca\e)±0.13%{stat.), 

(scale) ±0.70%(stat.), 
where we can see a clear reduction of the scale de- 
pendence at NLO. The statistical error coming from 
the numerical integration is negligible compared to 
the residual scale uncertainty. 

In the distributions the reduction of the theo- 
retical uncertainties also becomes apparent. The 
shaded areas in the distributions denote the change 
of the distribution when varying the scales Mw < 
fJ'R = tiF < 4: ■ Mw. All the distributions also 
include estimates of the pdf uncertainties. To esti- 
mate the pdf uncertainty we calculate the variables 
AX+jj^ and AX^^^ which can be derived from the 
set of eigenvectors that are included in the pdf set. 
The exact definition of these variables and a de- 
tailed discussion can be found in 72] . For this cal- 
culation we used the eigenvectors of the 90 percent 



aLo[/&] = 39.57 ■ 
aNLo[/&] =44.5r 



40 




200 
\i [GeV] 



Figure 1: Dependence of the LO and NLO cross sections 
pp W+{^ Vee+)W~{-^ fJ-~i^ii)jj at = 7TeV on 
renormalisation and factorisation scales. We use /i = fin = 
fip and vary between < M 4 Myir . 



confidence level set. Each phase space point enters 
the plots three times. Once with the central value 
of the pdf set and in addition for the two variables 
^^max ^i^d AX^^^ . However for the actual phase 
space integration, i.e. in the calculation of the to- 
tal cross section and the adaption of the integration 
grid during runtime, only the central value of the 
pdf set enters. In this way we are able to include 
an estimation of the pdf error into the histograms. 

In Fig. [2] the p± distributions of the two tagging 
jets are shown. Note that for the leading jet, the 
NLO corrections lead to an increase of the cross 
section in the first bins and a decrease in the high 
P± region. This can be explained by the possibility 
of radiating a third jet at NLO, which reduces the 
total available energy for the two leading jets and 
hence softens the leading jet energy. 

A similar behaviour can be observed for the lep- 
tons. In Fig. [3] we plot the distribution of the trans- 
verse momentum of the electron, and one can ob- 
serve that the NLO corrections slightly shift the 
curve towards smaller p± values. 

Fig. m shows the distribution of the azimuthal 
opening angle 0e+/j- between the two leptons. 
This distribution, relevant in the context of Higgs 
searches, peaks at 0e+^- = tt, in contrast to the 
case of leptons coming from a H W^W^ decay, 
where the largest contributions stem from (jj^+^i- 

HH. 
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Figure 2: Transverse momenta of the two tagging jets at 
yfs = 7TeV. The jets are ordered according to their p^. 
The bands describe the pdf and statistical error as well as 
the scale variations between A/^v 5: M 4 My^ ■ 



3.4- Relevance of the different virtual contributions 
Now we briefly discuss the impact of tlie virtual 
contributions of type A, B, and C described in Sec- 
tion 2. The numerical relevance of parts B and C is 
small: at the scale /i — 2-M\y, the absolute value of 
part B amounts to about 1% of the size of part A 
and part C to about 3%. The contributions com- 
ing from parts B and C have been systematically 
included in all the results presented in this paper. 

It is interesting to see whether the three types 
of virtual corrections differ in their kinematics. In 
Fig. Owe plot the individual contribution of Part A, 
B, and C to the p± distribution of the leading jet 
(// = 2 • Mw)- The three parts are all normalised to 
the next-to-leading order cross section at the scale 
/i = 2 • Mw in order to project out only their kine- 
matical behaviour. Parts A and B show a very sim- 
ilar shape whereas part C peaks at higher p± values 
and has a broader spectrum. As the three different 
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Figure 3: Transverse momentum of the electron, in pp 







at = 7 TeV. 
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Figure 4: Distribution for the opening angle between the 
leptons, 



parts do not necessarily have the same behaviour 
under a variation of the scales, the differences in 
shape will be different at a different scale. How- 
ever, as the overall sizes of part B and part C are 
small, we can infer that their contribution does not 
affect much the shape of the distributions shown in 
Section 13.31 

4. Conclusions 

We have calculated the NLO QCD corrections 
to the production of a W~^W~ pair in associ- 
ation with two jets at the LHC, including lep- 
tonic decays of the W-bosons with full spin cor- 
relations. The calculation was performed using the 
GoSam [HI] package for the virtual contributions. 
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Figure 5: Contribution of the three types of virtual correc- 
tions to the leading jet p± distribution. Each contribution 
is normalised to the next-to-leading order cross section. For 
a better visibility the distributions for part B and C have 
been multiplied with —10 and 10 respectively. For the scale 
we chose here fi = 2M-\^. See text for further details. 
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Appendix A. Numerical results at ampli- 
tude level 



in combination with MadGraph/MadEvent [53|, l56j 
and MadDipole [s^, [EHl for the real radiation and 
the phase space integration. The virtual amplitude 
also includes contributions with vector bosons at- 
tached to closed fermion loops (called part B) and 
contributions from quarks of the third generation in 
the loops (called part C), which have not been con- 
sidered in the previous literature. GoSam provides 
a setup in which these contributions can be easily 
filtered at the level of Feynman diagrams. We find 
that, at the scale = 2Mw, part B amounts to 
about 1%, while part C makes up about 3% of the 
size of the virtual contributions. The latter part 
also shows a different shape in its distribution, 
but as it is numerically small, the impact on the 
full distribution is not significant. 

Overall, we find NLO corrections of the order of 
10% for standard cuts at the LHC. The scale de- 
pendence is considerably reduced by the NLO cor- 
rections, from about 30% at LO to about 6% at 
NLO. Further phenomenological studies, including 
also vector boson fusion cuts, are in progress. 

Our results show that the package GoSam for au- 
tomated one-loop calculations, in combination with 
a suitable program for the real radiation part, is 
able to deal with non-trivial multi-leg calculations 
efficiently. 



To facilitate comparisons, we use the same phase 
space point as in [26|] but with the conventions 



Pi +P2 



i=3 



from 



Numerical values for the momenta pi . . 
Ref. [ii] are reported in Table [Ol 

The renormalization scale is fiji — 150 GeV, and 
for the points given in this appendix we use Mw = 
80.419 GeV, Tw = 2.141 GeV, Tz = 2.49 GeV, 
sinOw — 0.47138095. The numerical results for the 
ratios of the unrenormalised virtual amplitudes to 
the tree-level amplitudes are given in Table IA.2I 
The eight partonic configurations correspond to the 
basic partonic processes, listed in Section 2.2, from 
which all the others can be obtained via crossings 
and/or change of quark generation. 

For a more detailed comparison, we also provide 
the results for the finite parts of the unrenormalised 
amplitudes in the form a + bNf. The results, for all 
subprocesses, are contained in Table 1X^31 

Tables IA.4I and IA.5I contain two classes of dia- 
grams that have not been considered before, namely 
contributions in which vector bosons are attached 
to a closed fermion loops (Table E4|, and contribu- 
tions coming from the inclusion of the third quark 
generation in the loops (Table rATSj) . Only the par- 
tons accompanying the W^-boson pair are listed. 
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Px 


Py 


Pz 


pi 


500 


500 








P2 


500 


-500 








Pq/g 


214.48887016141776 


-27.06079802177751 


-98.519808378615 


188.59224795994947 


Pq/a 


54.23140701179994 


-31.13301620817981 


-7.9279665679114 


43.69128236111634 


Pl^e 


85.5312248384887 


-8.22193223977868 


36.16378376820329 


-77.0725048002413 


Pe + 


181.42881161004266 


-57.85998294819373 


-171.863734086635 


-5.611858984813 


Pm- 


82.84930107743558 


-65.90954762358915 


-49.89521571962871 


5.51413360058664 




381.47038530081545 


190.18527704151887 


292.042940984587 


-155.113300136598 



Table A.l: Kinematic point used for all the tables in this Appendix 



Part A 










a + hNf 


a 


b 


u, u, c, c 


-5.333333 


7.587051 


5.395241 




u, u, c, c 


14.575491 


-2.295061 


u, u, d, d 


-5.333333 


7.679374 


2.155533 




u, u, d, d 


10.222774 


-2.016810 


d, d, s, s 


-5.333333 


7.587051 


7.040878 




d, d, s, s 


16.296114 


-2.313809 


u, u, u, u 


-5.333333 


6.828238 


-2.016698 




u, u, u, u 


5.169238 


-1.796484 


d, d, d, d 


-5.333333 


6.733070 


3.214751 




d, d, d, d 


9.441580 


-1.556707 


d, u, c, s 


-5.333333 


7.587051 


-15.91575 




d, it, c, s 


-19.781800 


0.966511 




-8.666666 


-2.786885 


-4.712831 




u,u,g,g 


-4.660510 


-0.013080 


d, d,g,g 


-8.666666 


-2.836720 


-0.7309734 




d,d,g,g 


-0.678101 


-0.013218 



Table A. 2: Numerical results for the ratio of virtual over tree- 
level squared amplitudes summed over helicities and colors, 
A'^ = 4, in the same approximation as in Ref. [2^ . i.e. with- 
out contributions from 3rd generation quarks in the loops 
and without contrubutions where vector bosons are attached 
to closed fermion loops. Only the partons accompanying the 
W-boson pair are listed. 
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Part B 


e-' 


e-' 




u, it, c, c 








0.6750172 


u, u, d, d 
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d, d, s, s 
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u, u, u, u 








1.226365 


d, d, d, d 
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d, u, c, s 



















0.6780958 


d, d,g,g 
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